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Are there stable long-range ordered Fei x Cr x compounds? 
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The heat of formation of Fe-Cr alloys undergoes an anomalous change of sign at small Cr concentrations. This obser- 
vation raises the question whether there are intermetallic phases present in this composition range. Here we report the 
discovery of several long-range ordered structures that represent ground state phases at zero Kelvin. In particular we 
have identified a structure at 3.7% Cr with an embedding energy which is 49 meV/Cr atom below the solid solution. 
This implies there is an effective long-range attractive interaction between Cr atoms. We propose that the structures 
found in this study complete the low temperature-low Cr region of the phase diagram. 



Ferritic iron-chromium alloys are important materials for 
structural components in fusion and fission reactors. The need 
to understand their behavior under irradiation has recently 
motivated a significant amount of basic research the results 
of which have changed our understanding of these systems. 
One of the most surprising and challenging observations is the 
coexistence of short-range order and phase segregation in the 
same alloy: At low Cr-concentrations both experiment^— and 
simulations 8,9 have found a short range ordering tendency: Cr 
atoms seek to maximize the number of their Fe neighbors. 
At larger Cr concentrations, however, phase segregation is 
observed, i.e. Cr atoms favor Cr neighbors. Some insight 
into the origin of this peculiar behavior has been provided by 
several first-principles calculations 1 ^* which found a negative 
heat of formation at low Cr concentrations and an inversion 
of the heat of formation above approximately 10% Cr. These 
observations raise two important questions: (?) Are there inter- 
metallic phases in the Fe-Cr system which are missing in the 
phase diagram as we know it, i.e. is there long-range order? 
(if) Is the Cr-Cr interaction in the Fe host crystal purely repul- 
sive, i.e. do Cr atoms in Fe behave like a lattice gas, or is there 
an effective attraction between the Cr atoms? In the follow- 
ing we will demonstrate by means of density-functional the- 
ory (DFT) calculations that (i) there are several intermetallic 
phases at small Cr concentrations and that (if) some of these 
phases correspond to an effective long-range attractive inter- 
action between the Cr atoms. 

In order to motivate our approach we show in Fig. Q] a 
schematic plot which illustrates the consequences of both 
purely repulsive as well as attractive effective interactions be- 
tween the Cr atoms. For the following qualitative discussion 
we treat the Fe matrix as an effective medium into which the 
Cr atoms are embedded. The energy gained if a single Cr 
atom is removed from a bulk Cr crystal and inserted into a 
pure Fe matrix defines the heat of solution for a single (iso- 
lated) impurity, i.e. the embedding energy at infinite dilution 
which is indicated by the horizontal dashed lines in Figs. Q] 
and EJa). In the case of purely repulsive Cr-Cr interactions 
adding any other number of Cr atoms can only render the em- 
bedding energy per Cr atom less negative which leads to the 
shaded region in Fig. Q] In contrast, if attractive forces are 



present there will be configurations in which the embedding 
energy per Cr atom is more negative than the dilute limit as 
indicated by the solid circles in Fig. Q] In either case, since 
the alloy segregates for larger amounts of Cr, the embedding 
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FIG. 1. (Color online) Purely repulsive vs weakly attractive inter- 
actions. If the Cr-Cr interaction is purely repulsive the embedding 
energy per Cr atom for any possible configuration can only be less 
negative than the embedding energy at infinite dilution. If there is, 
however, some effective attraction between the Cr atoms, intermetal- 
lic phases exist and there are structures for which the embedding 
energy per Cr atom lies below the infinite dilution limit. 
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FIG. 2. (Color online) Calculated formation energies plotted (a) as 
embedding energy per Cr atom and (b) formation energy per atom. 
In (a) the structure with the lowest formation energy per Cr atom is 
marked by SI. In this structure the Cr atoms occupy a bcc sublattice 
three times as large as the underlying Fe bcc lattice. In (b) the most 
stable intermetallic phases SI and S2 are connected by thick solid 
lines. For compositions which fall on these lines, the system decom- 
poses into a two-phase mixture. For concentrations above 6.25% Cr 
the thick solid line connects S2 with pure Cr. 



energy must eventually become positive as more and more Cr 
atoms are added. Finally, it should be noted that in the case 
of purely repulsive interactions intermetallic phases are possi- 
ble but need not exist. However, if the interaction is attractive 
intermetallic phases must exist. 

In alloy thermodynamics, instead of considering the em- 
bedding energy, it is more common to describe a system in 
terms of the formation energy (identical to the mixing energy 
in the case of a random system), which is the energy per atom 
of the alloy compared to the fully separated components (both 
Fe and Cr). If one adopts this view the embedding energy 
at infinite dilution determines the slope of the dashed line 
in Fig. |2{b). If the Cr-Cr interaction is purely repulsive, all 
configurations with finite Cr concentration lie above this line, 
whereas at least some phases must lie below this line if there 
is some attraction. 

In several previous studie s 4,8,10 ' 11 the formation energies 
of a number of Fe-Cr configurations were calculated us- 
ing supercells derived from the conventional two-atom body- 
centered-cubic (bcc) cell (also compare Fig. |2). On the basis 
of these data one would have to conclude that the interaction 



between Cr atoms is purely repulsive. 

In the present study we have followed an alternative ap- 
proach. We systematically generated a large number of struc- 
tures in which the Cr atoms selectively occupy certain neigh- 
bor shells. No restrictions were imposed with regard to the 
symmetry of the Cr sublattice. Specifically we did not con- 
struct our supercells on the basis of the conventional two-atom 
body-centered-cubic (bcc) unit cell but considered multiples 
of the primitive cell. Any supercell employed in the present 
work contains exactly one Cr atom and is uniquely determined 
by an integer matrix, A, which relates the matrix of primitive 
cell vectors, h, to the vectors spanning the supercell, H, 



H = Ah with h 




(1) 



where the sum convention applies. Note that the determinant 
of A equals the number of atoms in the supercell. Using ap- 
propriate matrices A one can construct supercells with a vari- 
ety of different symmetries. Specifically supercells with face- 
centered cubic (fee) and simple cubic (sc) symmetries are ob- 
tained if one uses integer multiples of the matrices 



Ifcc 




and A*. 




(2) 



We generated a large set of different supercells focusing on 
structures with a depletion of Cr atoms in the first five neigh- 
bor shells. We then employed DFT to calculate the total ener- 
gies for these structures 19 and obtained the formation energies 
from 



AE = S FeCl - n Fe /(n F e + l)E Fe - E Cr 



(3) 



where np e denotes the number of Fe atoms in the supercell, 
EpeCr is the total energy for the supercell, £pe is the total en- 
ergy of the structurally identical pure Fe cell, and Ec r is the 
energy per atom obtained from a primitive anti-ferromagnetic 
Cr cell. If AE is divided by the number of Cr atoms in the 
cell one obtains the embedding energy per Cr atom, AE em b, 
whereas division by the total number of atoms yields the for- 
mation energy, AEf. Special care was taken to minimize the 
error in our calculations by varying the cell shape and volume 
to reduce strain effects. Formation energies were calculated 
exclusively as differences between fully equivalent supercells 
using fully equivalent fc-point grids. 

The key results of our calculations are compiled in Fig. [2] 
The embedding energy at infinite dilution has been obtained 
using a 108-atom fee supercell (0.9% Cr). For this struc- 
ture the neighbor shells of any Cr atom are depleted of other 
Cr atoms except for the 25th shell for which 12 out of 36 
neighbors are Cr. This corresponds to a Cr-Cr separation of 
4.24 an along the (411) direction where ao is the bcc lattice 
parameter. The embedding energy obtained for this cell is 
— 157 meV/Cr atom and henceforth is considered as the in- 
finite dilution limit. For comparison, the lowest concentra- 
tion of 1.9% considered in previous studies corresponds to a 
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structure in which the closest Cr-Cr separation is 3an. along 
the (300) direction and which yields an embedding energy of 
— 1 53 me V / Cr atom. 

Our most important discovery is that there are several struc- 
tures the energies of which lie below the line given by the in- 
finite dilution limit. According to the discussion above the 
location of these data points implies that there is a net at- 
tractive interaction between Cr atoms. The structure which 
shows by far the lowest energy per Cr atom is marked in Fig. [2] 
by the leftmost black circle (SI). In this 27-atom supercell 
with 3.7% Cr the Cr atoms are arranged on a bcc superlat- 
tice three times as large as the underlying bcc Fe lattice. The 
first nine neighbor shells are exclusively occupied by Fe atoms 
whereas the 10th shell is filled to one quarter by Cr atoms. 
The Cr-Cr pairs are aligned along the (111) direction and 
separated by 2.60 ar> The embedding energy for this struc- 
ture is —206 meV/Cr atom which is 49 meV/Cr atom below 
the infinite dilution limit. There are several other structures at 
concentrations smaller than 3.7% Cr with energies below the 
infinite dilution limit. They resemble the lowest energy struc- 
ture in that they display a depletion of Cr atoms at short range 
and partly or fully Cr filled shells at longer range. 

For concentrations larger than 3.7% Cr the structure with 
the lowest energy (S2 in Fig.|2| is based on a 16-atom base- 
centered monoclinic unit cell with lattice vectors 
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and 




(4) 



The S2 structure has an embedding energy of 
— 145 meV/Cr atom equivalent to a formation energy of 
— 9.1 meV/atom. In order to describe this structure on the 
basis of the conventional two-atom unit cell, one would have 
to employ a 16 x 16 x 8 (4096-atom) supercell. The S2 
structure described here is different from the intermetal- 
lic structure described in Refs. [l(] and [Tl] although both 
structures correspond to a Cr concentration of 6.25%. The 
structure described in Refs. [l0| and yjj, however, can be 
represented in a 2 x 2 x 2 (16-atom) cubic supercell in which 
the Cr atoms are arranged on a simple cubic lattice with a 
lattice constant which is twice as large as the lattice constant 
of the underlying bcc Fe matrix. Furthermore, its formation 
energy is only — 6.5meV/atom (compare the open triangle 
in Fig. |2] and Ref.O) and thus not as negative as the value 
calculated of — 9 meV / atom for the S2 structure. 

In order to determine the zero temperature phase diagram 
we consider the energy per the total number of atoms (for- 
mation energy) as shown in Fig. HJb). The infinite dilution 
limit is indicated by the dashed line and the stable phases are 
connected by thick solid lines, along which the system decom- 
poses into two-phase mixtures. The two relevant intermetallic 
ground state phases occur at 3.7% (SI) and 6.25% (S2) corre- 
sponding to the two structures described above. 

It is important to address the magnitude of the observed ef- 
fect. Even for the most optimal structures the energy gain per 
atom amounts to merely lOmeV/atom which corresponds to 
approximately 110 K. At higher temperatures the configura- 
tional entropy favors the solid solution. Due to the reduced 



atomic mobility at these low temperatures direct experimental 
observation of the ordering may be challenging. Neverthe- 
less our discovery suggests that the low temperature phase 
diagram of Fe-Cr alloys is more complex than previously 
thought. 

In summary, in this letter we have demonstrated that there 
are intermetallic ground state phases at the Fe-rich end of 
the Fe-Cr phase diagram with embedding energies as low 
as — 206 meV/Cr atom and formation energies as low as 
— 9meV/atom. The observation of embedding energies be- 
low the infinite dilution limit implies that there are effective 
long-range attractive interactions between the Cr atoms which 
are the result of many-body effects. In fact, attempts to cap- 
ture the observed behavior in terms of simple pair interaction 
models were unsuccessful. Extensive further study is required 
to elucidate the physical origin of this phenomenon. 
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partment of Energy by Lawrence Livermore National Lab- 
oratory under Contract DE-AC52-07NA27344 with support 
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gram. Generous grants of computer time through the National 
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Berkeley National Laboratory are gratefully acknowledged. 
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